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Abstract 

We investigate the local quantum uncertainty (LQU) between a block of L qubits and one single qubit 
in a composite system of n qubits driven through a quantum phase transition (QPT). A first-order QPT 
is analytically considered through a Hamiltonian implementation of the quantum search. In the case of 
second-order QPTs, we consider the transverse-field Ising chain via a numerical analysis through density 
matrix renormalization group. For both cases, we compute the LQU for hnite-sizes as a function of L and 
of the coupling parameter, analyzing its pronounced behavior at the QPT. 
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1. Introduction 

The interplay between quantum information the¬ 
ory and statistical mechanics has brought em erg ing 
connections between these research fields P, 0 0 . 
In particular, it has provided a deeper understand¬ 
ing about the role played by correlations in quan¬ 
tum phase transitions (QPTs). A seminal result 
in this direction is a link between the scaling of 
pairwise entanglement and QPTs in quantum spin 
chains ii- This has been further developed by 
introducing a distinction between the characteriza¬ 
tion of first-order and continuous QPTs 00]- For 
a block analysis, entanglement entropy has been 
found to be related to the central charge of the Vi- 
rasoro algebra associated with the conformal field 
theory behind the critical model More 

generally, it has been shown that quantum cor¬ 
relation measures such as provided by the quan¬ 
tum disco rd lllll are also able to identify quantum 
criticality [l2l [l3| . Remarkably, pairwise quantum 
discord may exhibit a more robust characteriza¬ 
tion of QPTs than pairwise entanglement in cer¬ 
tain cases. For instance, pairwise quantum discord 
between distant sites in a quantum chain may indi¬ 
cate a quantum critical point, while entanglement 
is absent already for very short distances IJ, [15 1. 


In addition, for finite temperatures, pairwise quan¬ 
tum discord is able to reveal the QPT by non- 
analyticities in its derivatives, while the pronounced 
behavior in two-qubit entanglement disappears for 
even small temperatures [l^ . 

In this work, we aim at investigating the behav¬ 
ior of the local quantum uncertainty (LQU) [I^ at 
quantum criticality. The LQU has been introduced 
as a quantum discord-like measure, which is primar¬ 
ily related with the skew information [3i 13 • 
particular, it plays a role in the characterization of 
quantum metrology protocols [I113- The behav¬ 
ior of LQU between pairs of spins in a quantum 
spin chain has been recently considered [2l|, [22 1. 
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Here, we generalize this previous analysis for sys¬ 
tems of dimension 2 x 2 by considering the LQU 
for blocks of arbitrary dimension D x 2 and also 
by discussing its finite-size behavior in both first- 
order and second-order QPTs. More specifically, we 
will evaluate the LQU between a block of L quan¬ 
tum bits (qubits) and one single qubit in a com¬ 
posite system of n qubits. For a first-order QPT, 
we will consider a Hamiltonian implementation of 
the quantum search, which is designed to find out 
a marked element in an unstructured search space 
of iV = 2" elements. By analytical evaluation, we 
will show that the LQU exponentially saturates to 
a constant value at the critical point as we increase 
the block length L. This saturation is found to 
be enhanced by the system size n. On the other 
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hand, at non-critical points, the LQU will be shown 
to vanish for large n. In the case of second-order 
QPTs, we consider the transverse-field Ising model 
with open boundary conditions. By implementing 
a numerical analysis via density matrix renormal¬ 
ization group (DMRG), we will show that the con¬ 
cavity of the LQU as a function of the block size 
L characterizes the QPT. For both first-order and 
second-order QPTs, we also consider the LQU as a 
function of the coupling parameter, showing that 
the LQU exhibits a pronounced behavior at the 
quantum critical point independently of the block 
sizes of L qubits. In particular, this pronounced be¬ 
havior is sensitive to n, showing a scaling behavior 
as we increase the size of the system. 


2. Local quantum uncertainty 


The uncertainty of an observable K in a quan¬ 
tum state p is usually quantified by the variance 
V{p,K) = TrpK^ — (TrpK)^. It may exhibit con¬ 
tributions from both quantum and classical sources. 
Quantum uncertainty comes from the noncommu¬ 
tativity between K and p, being quantified by the 
skew (not commuting) information [30 

I{p, K) = TrpK^ - Trp^/^Kp^/'^K. (1) 


Indeed, suppose p and K commute. Then, p and 
K have a common basis of eigenstates {|fc)}, which 
means that the uncertainty of K in an individual 
eigenstate \k) vanishes. Hence, a nonvanishing un¬ 
certainty V {p, K) is only possible if p is a classical 
mixing of {|fc)}. Therefore, the commutation of p 
and K implies that V {p, K) has a classical origin. 

The quantum uncertainty is intrinsically con¬ 
nected with the concept of quantum correlation. 
For example, let us consider a Bell state of two 
qubits, namely, \'tp) = (|00) -I- |ll))/-\/2, where 
{|0), |1)} denotes the computational basis. This is 
an eigenstate of the global observable cr^ 0 so 
there is no uncertainty on the result of a measure¬ 
ment of such an observable. On the other hand, 
the measurement of local spin observables is intrin¬ 
sically uncertain for the density operator \il}){il)\^ 
since an entangled state cannot be an eigenstate 
of a local observable. In particular, the variance 
V (p, K) for a local observable K will vanish if and 
only if the state is uncorrelated. 

The concept of quantum uncertainty can be ex¬ 
tended to mixed states. In this case, the skew in¬ 
formation I(p,K) vanishes if and only if p is not 


disturbed by the measurement oi K. If Lf is a 
local observable, the states left invariant by local 
measurement are the states with zero quantum dis¬ 
cord with respect to that local subsystem [3| ■ The 
quantum uncertainty on local observables is then 
intimately related to the notion of quantum dis¬ 
cord and, as shown in Ref. 0 , it can be used as a 
discord-like quantifier. We are now ready to define 
the local quantum uncertainty (LQU). Let p = pab 
be the state of a bipartite system, and let de¬ 
note a local observable on B {K is represented by a 
Hermitian operator on B with nondegenerate spec¬ 
trum F). The LQU as defined in [4], is given by 


Q{p) =minKrI{p,K^). 


( 2 ) 


Notice that Q is the minimum quantum uncertainty 
associated to a single measurement on subsystem 
B. If there is a if for which Q = 0 then there is no 
quantum correlation between the two parts of the 
state p. As proved in Ref. 17| , the LQU satisfies all 
the good properties of a discord-like measure. An 
analytical expression for Q can be obtained if we 
consider a bipartite D x 2 system. In this case 


Q{pab) = 1 ^max (Wab), (3) 


where Xmax is the maximum eigenvalue of the 3x3 
symmetric matrix W whose elements are given by 

(.WAB)ij = Tr[p]lg{lA 0 criB)PAB{lA 0 ^js)]- (4) 

In this work, we will consider a set of n qubits 
aligned in a chain, with the bipartition in subsys¬ 
tems A and B chosen as shown in Fig. [TJ 


n 


L 1 

Figure 1: Bipartition used to defined the subsystems A and 
B for the LQU evaluation. The size L of the block A is 
arbitrarily chosen and subsystem B is taken as one qubit. 


3. LQU for the quantum search 

The aim of the search problem is to find out a 
marked element in an unstructured list of N candi¬ 
dates. In a quantum setting, it is possible to solve 
the search problem with scaling \/]V, as proved by 
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Grover [^ . Here, we consider a Hamiltonian im¬ 
plementation through a quantum system composed 
of n qubits, whose Hilbert space has dimension 
= 2". We denote the computational basis by 
the set {|i)} (0 < J — 1). Without loss of gen¬ 
erality, we can assume an oracular model such that 
the marked element is the state |0). So the imple¬ 
mentation of the quantum search can be achieved 
through the projective Hamiltonian 

H{s) = (1 - s)(lL - |V'o)(^/’o|) + s(lL - |0)(0|), (5) 

where |'0o) = (1/V^) I*)j s denotes the 

normalized time 0 < s < 1. By preparing the sys¬ 
tem in its ground state at time t = 0 and by con¬ 
sidering an adiabatic dynamics, it evolves to the 
corresponding instantaneous ground state at later 
times. In particular, the system exhibits a first- 
order QPT at s = 1/2. The ground state energy in 
terms of the normalized time s reads 


Eis) = 


I - s)N 


( 6 ) 


with N = 1 — 1/N. For the ground state vector 
|'0(s)), we obtain 


IV'(s)) = V^|0) -b 


N-1 


( 7 ) 


where we have defined the quantities a(s) = 

i-i-(Af-i)fe= ’ “ i-i-(Af-i)fe21 = 1 — 


E{s 


Note that, in the thermodynamic limit n —>■ oo, the 
structure of the Hamiltonian implies that the LQU 
can only be non-vanishing at the quantum critical 
point, even though its scaling is nontrivial at finite 
sizes. This can observed from Eq. where both 
|0) and \tpo) are product states that become orthog¬ 
onal for n —>■ oo. In this limit, the ground state is 
|^/>o) for 0 < s < 1/2, with energy E{s) = s, while 
the ground state is |0) for 1/2 < s < 1, with energy 
E{s) = 1 — s. At s = 1/2 the ground state is de¬ 
generate. From Eq. d?]), |V'(l/2)) will be an equal 
superposition of |0) and |^/o) for n —>■ oo. It then 
follows that <5 = 0 everywhere except at s = 1/2. 

In order to determine the scaling at finite size 
n, we consider the density matrix p = \'ip{s))('4’{s)\ 
describing the system in the ground state, which 
can be written as 


pis) = 


abb 
b c c 


( 8 ) 


where b = y/a{s)c{s). As we trace out n' qubits 
of the system, the resulting partial density matrix 
p'{s) will be given by 


P'is) 


a' b' b' ... b'' 
b' d c' ... c' 


b' 


( 9 ) 


where a' = a -b (2"' - l)c, b' = b + (2"' - l)c, and 
c' = 2" c. We observe that p' is an N' x N' matrix, 
with N' = . Taking the square root from 

Eq. ([5]), we obtain 


V^ = 


Or 

bj. 


br 


br 

Cj' 


C’P 


br 

C-p 


Cp 



( 10 ) 


where we have defined ar = + X , b^ = 

-b X~[3~, and Cr = A+(/3+)^ -b X~{P~)'^, with 
r = ^4{b')‘^{N' - 1) -b (o' - c'{N' - 1))2, /3± = 

and A± = yJd±i^^/[l + [N'- 
l)(/3^)^). By rewriting Eq. (IlOp in block form, we 


get 


V^ = 


'A B B ... B 

B^ C C ... C 

B^ C C ... C 


( 11 ) 


where A, B and C are 2x2 matrices defined as 



( 12 ) 


Now we define Mi = yfd • (1® cr^), where Oi are the 
2x2 Pauli matrices. We now compute the matrix 
elements {WAB)ij = Tr[MiMj]. This can be ana¬ 
lytically performed for arbitrary size N, yielding 

Wip = 2 (6^ -b Crttr + AqbrCr + C^r) I 

W2,2 = 2 (crOr — 6^) ) 

W3,3 = al- ‘2-bl -b cl, 

Wap = IFpa = 2 [brUr - brCr + qbl - qcl) , 

( 13 ) 

with q = {N'/2 — 1) and the remaining matrix el¬ 
ements of W vanishing. We notice that the matrix 
W would be left unchanged if we interchange blocks 



























A and B, taking block A as the first qubit from the 
left and B as the L remaining qubits. In order to 
obtain the LQU, we have to find out the largest 
eigenvalue of Wab and use it into Eq. Since 
four entries of the matrix W have vanished, this can 
be done analytically as well. 



Figure 2: (Color online) LQU between a block of length L 
and a single qubit at the quantum critical point s = 0.5. 
An exponential saturation as a function of L is observed, 
with the saturation value enhanced by the system size n. 
The inset shows the curve for n = 40, which approaches the 
asymptotic form and can be fit by Q(L) = 1/2 — ae~^^ , with 
a ^ 0.29 and b ss 0.71. 

By taking the LQU between a block of length L 
and a single qubit, we obtain an exponential sat¬ 
uration of the LQU at the quantum critical point. 
This is illustrated in Fig. [2j where we plot the LQU 
as a function of L for s = 1/2. As we can see, the 
saturation value is enhanced by n and L, attaining 
(5 = 1/2 for n —>■ oo and L —>■ oo. In this regime, 
the composite system AB is described by a pure 
state, with the LQU equivalent to bipartite entan¬ 
glement as measured by the linear entropy [ 13 . By 
computing the linear entropy S{pb) = 2 (l — Trp^) 
from Eq. © for n —>■ oo and L —>■ oo, we directly 
obtain S{pb) = 1/2, in agreement with the thermo¬ 
dynamic limit in Fig. [51 On the other hand, as dis¬ 
played in Fig.jSl the behavior of the LQU displays a 
” tilde” form if the system is driven to a non-critical 
point. In this case, as the total number of spins n 
grows, the curve tends to vanish independently of 
the block length L. 


The behavior of the LQU as a function of s for a 
fixed block size L is also remarkable. Since we have 



Figure 3: (Color online) LQU between a block of length L 
and a single qubit at the non-critical point s = 0.49. The 
LQU tends to vanish as we increase the size n of the system 
(see Inset). 

a first order QPT, the LQU itself shows a sharp 
behavior at s = 1/2, as shown in Fig. j?] 



Figure 4: (Color online) LQU between a pair of qubits (L = 
1) as a function of the normalized time s in the ground state. 
The system shows a sharp behavior at the critical point as 
a function of the system size n. In particular, exponential 
convergence of the LQU towards its thermodynamical value 
is obtained for each value of s. 


Indeed, this is exhibited for a pair of qubits for 
systems of distinct total number n of qubits. No¬ 
tice that, for L = 1, we obtain (5 = 1/3 in the limit 
n ^ oo. This value comes from the structure of 
|■^/(l/2)) in Eq. © for n —>■ oo, which is an equal 
superposition of the states |0) and |t/o)- In partic¬ 
ular, this maximum value of Q, which is obtained 
at the critical point s = 1/2, strongly depends on 
L. This is illustrated in Fig. [SI where it is shown 
that the maximum value of the LQU increases along 
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with the total number n of qubits in the system and 
L for the block, approaching Q = 1/2 in the limit 
of both L and n approaching infinity. 



Figure 5: (Color online) The maximum of the LQU for dis¬ 
tinct sizes L of blocks. 


4. The transverse field Ising model 

Let us consider now the Ising spin chain in a 
transverse magnetic field, whose Hamiltonian is 
given by 

n 

= (14) 

i=l 

with open boundary conditions assumed. Without 
loss of generality, we will set the energy scale such 
that J = 1. This Hamiltonian is Z 2 -symmetric and 
exhibits a second-order QPT from ferromagnetic to 
paramagnetic state at A = 1 [2^, [2^ . It can be ex¬ 
actly diagonalized by mapping it to a spinless free 
fermion model with single orbitals. However, in¬ 
stead of this solution, it is very convenient to con¬ 
sider a treatment based on density matrix renor¬ 
malization group (DMRG) [13, since (L-|-l)-point 
correlators and partial density matrices will be the 
objects of interest for the computation of the LQU 
for block systems. In the DMRG scenario adopted 
here, we truncate the system at each renormaliza¬ 
tion step of the density matrix, keeping a number 
M = 30 states for chain sizes up to 48 sites. This 
is well justified in the transverse-field Ising model, 
since the neglected states in the DMRG procedure 
have probabilities of the order of 10“^°. 



Figure 6: (Color online) LQU as a function of the block 
size L in the transverse-field Ising model at quantum critical 
point A = 1. 


We consider the behavior of the LQU as a func¬ 
tion of the block size L for A = 1, A = 0.8, and 
A = 1.2, which are exhibited in Figs. [51 [71 and [51 
respectively. As we can see. Fig. [51 shows that, at 
the critical point, the concavity of the LQU is un¬ 
defined. We can also observe that the LQU reaches 
its maximum value {Qmax ~ 0.41) independently 
of the number n of sites. This is indicated by the 
dashed line in Fig. [51 On the other hand, concavity 
is well-defined at off-critical points. Below the crit¬ 
ical point, LQU is a convex function of the block 
size L, as shown in Fig. [71 Oppositely, LQU is a 
concave function of L above the critical point, as 
displayed in Fig. [51 



Figure 7: (Color online) LQU as a function of the block size 
L in the transverse-field Ising model at A = 0.8 (below the 
critical point). 
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Figure 8: (Color online) LQU as a function of the block size 
L in the transverse-field Ising model at A = 1.2 (above the 
critical point). 

The quantum critical point can also be revealed 
as a function of the coupling parameter A for a fixed 
length L. By varying A, we can observe a sharp be¬ 
havior of the LQU around A = 1 for finite systems. 
This is illustrated in Fig. [HI where the LQU between 
a block of L = 10 sites and a single site is shown 
for distinct lengths n of chains. The inset in Fig. [ 5 ] 
shows the first derivative of the LQU with respect 
to parameter A. As a typical characterization of the 
QPT by a quantum discord-like measure [l^ . the 
critical point is identified by the first-derivative of 
the LQU with respect to A. 



Figure 9: (Color online) LQU as a function of the coupling 
parameter A in the transverse-field Ising model for distinct 
chain sizes n. Inset: First derivative of Q with respect to A, 
displaying the characterization of the the QPT. 

In Fig. [TUI 'we show the LQU and its first deriva¬ 
tive in the (A, L) plane. It is observed that around 


the critical point (A = 1), the value of the corre¬ 
sponding LQU on the color scale is approximately 
0.4. The first derivative of LQU is also shown in 
Fig. 1101 The sharp behavior of this quantity char¬ 
acterizes the QPT. The minimum value of the first 
derivative has a symmetry behavior as a function 
of block size L. 



Figure 10: (Color online) LQU and and its first derivative 
as a function of the coupling parameter A and the block size 

L. 


5. Conclusion 

In conclusion, we have investigated the scaling 
properties of the LQU at first-order and second- 
order QPTs. We have considered the behavior 
of the LQU in terms of either the block size L 
or the coupling parameter inducing the QPT. In 
both cases, the QPTs are precisely identified by 
the LQU. The scaling of the LQU as a function of 
the size of the block opens the possibility of search¬ 
ing for connections between quantum correlations 
and finite-size properties of critical systems. As a 
future step, it would be interesting to investigate 
possible universal properties of LQU, e.g. analyz¬ 
ing its behavior for quantum critical spins chains 
belonging to distinct universality classes. A rele¬ 
vant point would be a possible scaling related to 
the central charge of the Virasoro algebra behind 
the critical chain. Other points of interest would 
be the investigation of the quantum critical behav¬ 
ior for finite temperatures as well as the effect of 
boundary conditions over the scaling of the LQU 
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for systems D x 2. These investigations are left for 
further research. 
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